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SCIENTIFIC  BESDLTS  &ND  PBOGBESS  DURING  THE  PAST  12  HONTRS 


THEOBETICAL 

The  main  theoretical  results  of  our  past  year’s  work  is  con¬ 
tained  in  the  paper  with  T.Tsai,  ’Time  Dependent  Pulse  Am¬ 
plification  in  a  Three  level  Gas*.  The  abstract  is  given 
below. 


The  time  dependent  density  matrix  equations  have  been  sblved 
^  in  conjunction  with  the  field  amplitude  equations  for  a 
three  level  optically  pumped  gas.  The  solutions  are  valid 
in  the  limit  in  which  the  pumping  pulse  is  longer  than  the 
dephasing  time  of  the  system  and  true  coherent  effects  may 
be  neglected.  It  is  assumed  that  the  FIB  field  is  smaller 
than  the  pump  field,  and  saturation  effects  associated  with 
this  field  are  absent.  In  the  limit  in  which  the  pump  pulse 
is  resonant,  particularly  simple  results  for  the  FIR  field 
amplification  nay  be  obtained.  The  solutions  exhibit  the 
dependence  of  the  FIB  output  pulse  on  the  shape  of  the  input 
pulse.  Further,  as  a  consequence  of  non-linear  saturation 
effects  on  the  pump  pulse  transition,  the  FIB  pulse  is  ini- 
tialy  created  with  two  peaks,  one  associated  with  the  lead¬ 
ing  edge  and  one  with  the  trailing  edge  of  the  pump  pulse. 
As  amplification  on  the  FIB  transition  and  pump  absorption 
increase,  these  two  peaks  merge.f^ 


Enclosed  is  a  copy  of  the  paper  that  has  been  forwarded  to 
the  IEEE  Journal  of  Quantum  Electronics. 


EIPEBIHEHTAL 


Experimental  work  over  the  past  year  has  been  concerned  with 
the  detailed  measurement  of  breakdown  characteristics  of  ar¬ 
gon  gas.  It  had  been  hoped  that  b,y  measuring  the  time  de- 
pendennce  of  the  breakdown  in  a  gas,  that  certain  'finger¬ 
prints'  would  appear  that  independent  of  the  statistical 
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nature  of  laser  induced  gas  breakdown.  Several  interesting 
aspects  of  laser  induced  gas  breakdown  were  observed,  some 
of  which  were  not  reported  previously  in  the  literature,  and 
it  was  learned  that  the  only  characteristic  we  could  isolate 
that  did  not  vary  statistically,  was  the  slope  of  beam  at¬ 
tenuating  plasma  as  a  function  of  pressure.  It  had  been 
hoped  that  this  pressure  dependence  would  be  sensitive 
enough  for  this  technique  to  be  employed  as  a  diagnostic 
tool;  however,  this  dees  not  appear  to  be  the  case.  In  ad¬ 
dition  to  obtaining  the  first  laser  induced  gas  breakdown  in 
argon  in  which  the  temporal  behavior  of  the  breakdown  was 
monitored,  this  was  done  with  a  single  mode  laser  so  that  a 
more  accurate  value  of  threshold  could  be  obtained.  These 
results  are  being  prepared  for  publication  in  the  Journal  of 
Applied  Physics.  In  addition,  we  have  discovered  a  rela¬ 
tively  simple  method  to  measure  the  divergence  angle  of  a 
carbon  dioxide  laser  beam  using  burn  patterns  on  uncoated 
high  speed  polaroid  film.  This  depends  upon  the  existence 
of  a  gaussian  beam,  and  upon  the  distinct  two  level  sensi¬ 
tivity  of  the  film  to  the  burn  pattern.  He  hope  to  prepare 
these  useful  results  in  a  note. 

PKEVIOUS  ACCOHFLISHtlENTS 

Previous  accomplishments  have  been  associated  with  the  de- 
velopcment  of  kinetic  equations  for  polyatomic  gases,  and 
the  solution  of  various  boundary  value  problems  of  specific 


relevance  to  the  space  prograo.  (lore  recent  work  has  teen 
concerned  with  the  interaction  of  radiation  and  gases  in  a 
three  level  system,  and  equations  have  been  derived  from  the 
density  matrix  that  are  rate  equations  retaining  nonlinear 
power  broadening  and  multiphoton  effects.  These  equations 
have  been  solved  in  the  case  in  which  the  pumping  field  is 
stronger  than  the  FIP.  field,  and  such  solutions  should  be  of 
interest  to  experimentalists  concerned  with  the  optimization 
of  of  FIB  optically  pumped  lasers. 
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bv 

Tsunq-Minq  Tsai  and  T.F.  Morse 

Division  of  Enqineering,  Brown  Dniversitv, 
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ABSTRACT 


The  density  matrix  equations  have  been  solved  in  conjunc¬ 
tion  with  the  time  dependent  field  amplitude  equations  for  a 
three  level  optical Iv  pumped  qas.  The  solutions  are  valid 
in  the  limit  in  which  the  pumping  pulse  is  longer  than  the 
deohasing  time  of  the  system  and  true  coherence  effects  may 
be  neglected.  In  order  to  obtain  closed  form  solutions,  it 
is  assumed  that  the  FIR  field  is  smaller  than  the  pump  field 
and  saturation  effects  associated  with  this  field  are  ab¬ 
sent.  When  the  pump  pulse  is  resonant  with  the  absorbing 
FIR  transition,  particularly  simple  solutions  for  the  time 
dependent  FIR  field  amplification  are  obtained.  These  solu¬ 
tions  exhibit  the  specific  dependence  of  the  FIR  output 
pulse  on  the  shape  of  the  input  pulse  and  display  certain 
pulse  narrowing  features. 
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1.  INTRODUCTION 


There  has  been  much  recent  interest  In  both  experimental 
as  well  as  theoretical  studies  of  optically  pumtjed  qas  laser 
systems,  with  particular  emphasis  on  FIR  wayelengths. ^ 
excellent  reyiew  is  presented  in  reference  14.  Many  of 
these  systems  operate  most  efficiently  in  the  pulsed  mode 
and  our  analysis  will  be  restricted  to  a  fuller  understand¬ 
ing  of  the  nature  of  the  transient  response  of  such  systems, 
as  perhaps  best  typified  by  the  CH^F  system  pumped  with  a  CO2 
laser.  The  physical  model  we  use  is  that  of  the  three  leyel 
atom  or  molecule  that  has  been  extensiyely  inyestigat- 
ed. 15-18 

The  basic  framework  of  this  description  of  transient  be¬ 
havior  in  an  optically  pumped  gas,  is  that  of  the  density 
matrix  equations  in  the  limit  in  which  the  time  characteris¬ 
tic  of  the  pumping  pulse  is  longer  than  the  relaxation  times 
of  the  FIR  system.  In  essence,  we  are  examining  the  tran¬ 
sient  solution  to  the  equations  presented  in  references  2 
and  22.  Our  formulation  of  the  appropriate  equations  in  no 
way  introduces  any  new  physical  component;  however,  the  con¬ 
venient  form  in  which  thev  are  cast  permits  us  to  obtain  new 
solutions  for  the  time  dependent  behavior  of  optically 
pumped  systems.  In  that  which  follows,  we  shall  restrict 
our  attention  solely  to  homogeneously  broadened  systems. 

With  regard  to  earlier  solutions  obtained  for  the  tran¬ 
sient  behavior  of  optically  pumped  three  level  systems,  we 


note  that  orevious  authors  have  postu.late(1  phenomenological 
rate  equations  that  omit  certain  essential  aspects  of 
these  svstems. 5-9  features  often  neglected  in 

these  descriptions  are  associated  with  non-linear  power 
broadening  on  both  the  pump  and  the  FIR  fields,  and  the  mul¬ 
tiphoton  or  Raman  effects  that  may  be  present.  These  latter 
transitions  are  important  for  the  case  of  off  resonant  pump¬ 
ing  of  FIR  systems,  especially  since  pump  detuning  of,fers 
the  possibility  of  tunability  of  FIR  frequencies  over  a 
fairly  wide  range  of  wavelengths. 18-20  Solutions  to  simpli¬ 
fied  versions  of  three  level  rate  equations  that  neglect 
these  effects  can  yield  proper  results  onlv  in  certain  lim¬ 
its. 

In  the  following  section  we  briefly  develop  the  equations 
used,  and  in  section  3  we  present  our  general  solution.  In 
section  4,  certain  asymptotic  limits  are  discussed,  and  a 
comparison  with  some  experimental  results  is  made. 

2.  THREE  LEVEL  RATE  EQUATIONS 

2.1  NON-LINEAR  RATE  EQUATIONS 

From  the  density  matrix  equations  we  mav  properly  derive 
rate  equations  that  retain  important  power  broadening  terms 
and  multiphoton  effects  necessarv  for  a  description  of  the 
amplification  of  an  optically  driven  FIR  pulse. 1^“15  The  re¬ 
striction  in  this  derivation  is  that  the  pump  pulse  be  long- 


2 


er  than  the  dephasinq  time  of  the  system,  so  that  we  may  re- 
aard  all  of  the  off  diaqonal  terms  of  the  density  matrix  as 
synchronized  to  the  instantaneous  value  of  the  electromaq- 
netic  field. 22  Thus,  a  complete  rate  equation  description  is 
appropriate  in  all  situations  in  which  true  coherence  ef¬ 
fects  are  absent.  A  schematic  of  the  physical  system  we 
wish  to  analyze  is  shown  in  fiqure  1,  and  a  description  of 
the  detuninq  parameters  is  given  in  fiqure  2.  The  actual 

system  we  wish  to  describe  is  naturally  richer  in  detail 

schematically 

than  the  three  level  model  we  employ,  and  fiqure  SAillus- 
trates  some  of  the  true  complexities  of  a  typical  FIR  sys¬ 
tem.  These  problems  will  be  discussed  below  when  we  calcu¬ 
late  physical  parameters  to  make  comparison  with  some 
experimental  results.  The  consequences  of  averaqinq  over  M 
and  K  levels  have  been  quantitatively  discussed  in  a  steady 
state  analysis.'^ 

On  a  time  scale  longer  than  the  true  coherence  time  of 
the  particle  internal  states,  the  nature  of  the  interaction 
between  radiation  and  matter  may  properly  be  described  by  an 
appropriate  line  shape  function,  Q.  This  contains  informa¬ 
tion  on  the  response  of  the  system  to  resonant  and  near  re¬ 
sonant  radiation  and  may  incorporate  important  nonlinear 
broadeninq  effects  and  multiphoton  transitions.  In  this 
limit,  the  density  matrix  equations,  or  rather  the  equations 
for  particle  number  difference,  reduce  to  the  following.  We 
have  assumed  that  the  total  number  of  particles  in  states  1, 


2  and  3  is  constant 


3t 


32 


13 


Yi  Y2 

^3  ^4 


32 


32 


■  ^13 


2  2 
^^32  -  2^13 


"^^32 


'»^13 


ImQ 


32 


^13 


(1) 


Relaxation  coefficients  Yj^j  are  presented  in  Appendix  1, 

’^ii  “  ^ii  ”  ’^^ii  is  the  equilibrium  state  differ¬ 
ence  poDulation.  The  complex  line  shaoe  function  0  is  re¬ 
lated  to  the  susceptibility  in  the  followinq  manner. 
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The  Rabi  frequency  is  given  in  terms  of  the  electromagnetic 

field,  =  y  iiEi-j/fi,  and  standard  notation  of  reference 

1  is  employed.  The  equations  for  the  pump  and  FIR  field  am¬ 
plitudes  are 

(  Ir  *  c  i-l 

where  Y^  g^d  Ip  are  phenomenological  loss  terms.  Sponta¬ 
neous  emission  is  neglected.  It  may  be  incorporated  in  our 
solutions  as  an  initial  condition  on  the  FIR  field  intensi¬ 
ty.  The  equations  for  pump  (  ♦p)  and  FIR  phase  (  ^g)  are 
then  qiyen  by  the  following. 
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where  X  and  X  refer,  respectiyely ,  to  the  real  and  im¬ 
aginary  parts  of  the  complex  susceptibility.  These  expres¬ 
sions  are  algebraically  quite  Inyolyed  when  presents'!  in 


terms  of  the  line  shape  function  (Appendix  1)  and  equation 
2.  To  simplify  these  relations,  we  proceed  as  follows. 
Notinq  the  detuning  parameters  of  fiqure  2,  we  introduce  in 
Appendix  2  and  in  figure  4,  complex  detuning  parameters  that 
permit  the  collisional  dephasinq  frequencies  and  the 
FIR  and  pump  field  detuning  to  be  combined  as  complex  val¬ 
ues.  We  now  define  saturation  parameters  that  are  deoendent 
upon  detuning  as  well  as  upon  collisional  parameters. 
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With  these  quantities,  dimensioTiless  pump  and  FIR  field  in¬ 
tensities  may  be  introduced. 


S  =  B32^“i3  “i2  = 
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The  complex  line  shape  function  0  may  now  be  written  in  a 
more  compact  form  as  follows. 
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with  A  a  e^"^  +  el-‘*’i3  +  ei*>’32.  From  figure  4  we  see  that  <!• 
is  the  sum  of  the  three  detuning  angles,  ♦  ®  ®  ,  and 

0  =  ♦32  + 


-  i 

S  +  ■  *32^  P 

“12  ° 

A 

S  P  +  ■  '*’13^ 

0 

5 


We  now  are  able  to  write  the  population  difference  equa¬ 


tions  as  follows, 
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The  eouations  for  the  FIR  and  pump  intensities  and  phases 
are  given  as 
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Fiquations  9-11  must  be  solved  with  initial  conditions  speci¬ 
fied  on  the  state  densitv  differences  and  the  input  pump 
pulse  shaoe.  We  mav  take  the  initial  conditions  on  the  FIR 
field  to  be  either  background  black  body  radiation,  or  a 


weak  FIR  pulse  that  we  wish  to  amplifv.  These  equations  are 
still  formidable  and  nonlinear^  and  it  will  be  necessary  to 
consider  solutions  in  the  limit  in  which  saturations  effects 
on  the  FIR  pulse  mav  be  neqlected.  This  is  done  in  the  fol¬ 
lowing  section. 


2.2  RATE  EQUATIONS  FOR  S/P  «  1 

To  proceed  further  toward  our  goal  of  a  solution  with  a 
relative  degree  of  analytical  simplicity,  it  is  necessary  to 
make  one  additional  assumption  of  significance;  i.e.,  that 
the  ratio  of  the  dimensionless  FIR  field  to  the  dimension¬ 
less  Dump  field  is  small,  S/P  <<  1,  where 
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Tvpical  values  are  ^^3  «  .05  D,  and  ^^2  =  1  D,  so  that  if 
we  set  an  upper  limit  of  S/P  »  .1,  we  have  then  the  re¬ 
striction  that  Ijj/ip  <  2.5  X  10"^.  This  ratio  can  be  larger 

for  >  ®s*  Since  P  decreases  as  S  increases,  this  de¬ 

scription  will  rigorously  hold  only  for  certain  lengths  of 
amplifier.  For  lengths  such  that  saturation  in  the  FIR 
field  plays  a  role,  our  results.  may  be  only  qualita¬ 

tive.  With  this  additional  linearization  of  the  FIR  field. 


equations  9-11  become 
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tainina  power  broadening  only  on  the  pump  line  that  all 
eouations  (exceot  the  equation  for  the  FIR  field)  are  inde¬ 
pendent  of  S.  This  is  a  natural  consecnience  of  our  neglect 
of  FIR  saturation  effects.  We  have  also  made  a  standard 
transformation  to  a  retarded  time  variable,  t  «  t  -  z/c  and 
have  neglected  diffusion  effects  in  the  equation  for  state 
densitv  differences.  Equations  13-14  are  rate  eouations 
that  retain  nonlinear  power  broadening  effects  associated 
with  the  pump  pulse  as  well  as  contributions  from  multipho¬ 
ton  processes.  In  that  which  follows,  we  shall  obtain  solu¬ 
tions  to  equations  13-18  in  appropriate  limits. 


3.  SOLUTION  OF  THE  EOUATIONS 
Prior  to  considering  a  solution  for  the  pump  amplitude 
and  phase,  which  is  the  kev  to  the  solution  to  our  system  of 
equations,  we  introduce  the  following  complex  transfer  func¬ 
tions  whose  use  simplifies  the  FIR  and  pump  equations. 

•  J  1"  -  i  -  ♦p(o.T))  (15, 


A^Cz.t) 


A^Co.t)  =  i  In 


S(Z,T) 

S(o,t) 


-  i  [4>5(z,t)  -  ♦g(o,T)] 


(20) 


With  these  relations,  equations  15-18  may  be  witten  more  compactly 
as  follows. 


-Y 

=  -JB. 
2 


3z 


„  .  -ie 


P  +  e 


ie 


(21) 


P  ♦  e 


P 


(22) 


'’’he  three  absorption  cross  sections  Oj,  and  03,  are 

defined  in  Appendix  2.  In  order  to  integrate  ecruation  21, 
it  is  necessary  to  solve  for  r^j.  to  facilitate  this,  we 
assume  that  all  of  the  T-,  times  are  eoual,  and  that  all  of 
the  T2  times  are  equal.  Thus,  as  can  be  seen  from  Aooendix 
^11  “  ^22  *  ^33  =  l/^lf  ^2  “  ^4  =  0,  Y  =  Yj  =  Y3  1/Tj. 

Y  =  ■^^“32“l2  ♦l2''“l3^  “  "‘^'^2  (23) 


where  Tj  is  the  collisional  deohasinq  time. 

At  this  point,  it  is  convenient  to  define  4>  ,  a  pump 
Dulse  area  that  retains  the  effects  of  collisional  deactivi- 
ation. 


^^PC2,t)Ydt 


(24) 


If  we  now,  in  addition,  define  a  density  function  p 
that 


.  f  ' 

J  .OB 


-  t) 


rjjPYdt  .  I  x,jd» 
*0 


♦ 


then  equations  13  and  14  take  on  extremely  simple  forms 

^^32  -  -  o  Yd 

—  *  ^'32  ■  "''32  •  1  '■13'^ 


such 


(2S) 


(26) 


^  *  ^’-13  °  ”^3  ■  ’'’■13'’ 


(27) 


These  may  be  formally  inteqrated  to  yield. 


•■32  ■  '■32  * 


(28) 


-  10  - 


g  -  * 


Further 


2 


When  these  results  are  substituted  into  eouations  21  and  22, 
we  mav  integrate  to  obtain  the  following  exoressions. 


A„(z,t)  -  A^(o,t)  =  o,r°  e‘^*13 


1  13 


5,  [  *  Pe'*13  M 


d4> 


CIO) 


A  (z.T)  -  A  ro 

*  S 


^  I.  -  |^|L 

^  o  P  +  e^® 


d*' 


(31) 


Equation  30  and  31  are  formal  results.  Not  only  must  we  ob¬ 
tain  P  and  O*  /  3  2)^  as  functions  of  ♦  to  perform  the 
integration,  we  must  somehow  obtain  ♦  as  a  function  of  the 
pump  pulse  P,  or,  equivalently,  ?  =  ?(♦),  In  order  to  pro¬ 
ceed  further,  we  make  the  additional  assumption  that  the 
cavity  loss  parameter  for  the  pump  pulse,  may  be  neg¬ 

lected.  Multiplvinq  equation  16  by  e“^  ^  ^  and  integrat¬ 
ing,  we  find  that 


iil 

3zJt 


-OP 


(32) 


where  o  =  -2  cos  4)^^  This  permits  the  integration  of 
equati  n  30. 


„(z,t)  -  A^Co.t)  «zr°  o,e"^*13  ♦  (o,/o)e"^*13  (♦  -  ♦  J 


We  now  mav  seoarate  the  real  and  imaqinary  parts  of  the 
above  solution  to  obtain  the  pump  amplitude  and  phase  as 
functions  of  position  and  ♦  . 


In 


Pfz.Tl  ^“13  “/is"  +13 

- - - -  *  - - - - -  + 


P(o,t) 


(*o  -  ♦) 


(34) 


13 


(35) 

♦p(z.T)  -  4p(o.T)  =  -  CUi3(i>pr°3z/fiej:Yj2)  sin  cos  ^  (^^-Otan 

Equations  34  and  35  reoresent  a  solution  to  the  behavior 
the  pumo  pulse  if  ♦  =  4>CP)  is  known.  As  previously  indi¬ 
cated,  we  will  examine  this  relationship  in  the  limit  in 
which  the  time  characterizinq  the  pump  pulse  half  wi'^th  is 
longer  than  either  the  Tj  or  Tj  relaxation  times.  This  ap¬ 
proximation  corresponds  to  many  optically  pumped  FIR  laser 
systems  pumped  by  a  CO2  TEA  laser. 

For  the  case  in  which  ‘'pAi  >>  Ir  we  may  take  the  limit 
in  which  t  becomes  larqe,.  Thus,  the  exponential  in  equa¬ 
tion  24  mav  be  replaced  by  a  6  function,  and 

4T, 

♦  =  1  P  =  P  (36) 

From  equation  34,  we  obtain  the  followinq  result. 

P/Po  =  exp[7  X  (p^  .  p^ 


12  - 


This  simple  ecaiation  defines  the  spatial  and  temporal  varia¬ 
tion  of  the  pump  pulse  as  a  function  of  the  cross  section 
o,  the  initial  state  densities  and  the  ratio  of  re¬ 
laxation  times  or  eouivalentlv,  y  /  ir  .  The  depen¬ 

dence  of  P  upon  the  retarded  time  appears  onlv  throuqh 

^o^0»  )  f  the  initial  condition  on  the  pump  pulse.  Equation 

37  will  be  compared  with  some  exoerimental  results  below. 
The  solution  for  P(  ♦  )  to  be  substituted  into  our  exoression 
for  the  FIR  field  is  qiven  bv  equation  36.  The  remaininq 
ouantitv  to  be  determined  in  order  to  be  able  to  integrate 
equation  31  is  p(*K  We  note  the  singular  nature  of  the 
limit  in  which  Y  becomes  large  bv  writing  equation  27  in 
the  following  manner. 


^  (>-13  -  '•l3>  -  ’"'is'’ 


f  !!ul 

f  M.1 

[  3*  Jz 

1  3y  J 

Dif f eren tating  equation  24  with  respect  to  ^  ^  and  substi¬ 
tuting  into  the  above  relation,  we  find 

IP-r-  ■?  *1  -^''■13  •  '13’  '  ^^3’'  O’) 

In  the  limit  of  large  ^  ,  the  left  hand  side  of  equation  39 
vanishes,  and  we  obtain  an  algebraic  dependence  of  p  upon 
♦  .  ’’’his  merelv  indicates  the  limit  in  which,  for  the  Dump¬ 
ing  pulse  longer  than  the  relaxation  times  of  the  svstem, 
the  state  densities  respond  instantaneously  to  changes  in 
the  value  of  the  pump  field.  Thus, 


13  *13 


*  ♦) 


♦■"i/a  ♦  ♦) 
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We  now  have  the  relations  for  PM  ),  Pf  *  > ,  an<l  )  to 

be  substitute'!  into  equation  31.  The  integrations  may  nov; 
be  explicitly  carried  out. 


“32"s 


^“l3“p  *13 


P  e 

P  + 
o 


(41) 

Resolving  equation  41  into  real  and  imaqinarv  components,  we 
obtain  solutions  for  the  FIR  intensity  and  phase. 


where 


ffefi I 

U, 

13  p 

2 

^32*^5 

+  —2 -  (a  tan  41^2  "  4'l3) 

^13“p 


+  ij 
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1  *  ^0^1 


2 

♦  _(z.T)  -  45(0, T)  =  (a  *  1)  tan‘^  f - 

u  (1)  11  +  XqX/ 


2 

1  *^32***s 

j  -  (a  tan  *^2  '  *13^ 

u, 


fX?  +  I't 


»^13“p 


x^  +  ij 


(42) 


f43) 


„  P  +  cos  6 

^  “  siF"e  (44) 

-  J_  ^11 

*  “  2Y  032  cos  423 

These  solutions  have  required  the  following  assumptions. 
1.  'Phe  pumping  pulse  tp  is  longer  than  either  of  the  relax¬ 
ation  times  T^  or  T2.  2.  The  cavity  losses  on  the  pump 

pulse,  have  been  neglected.  3.  '’’he  total  population  of 

the  three  states  has  been  taken  as  constant.  4.  We  have  as- 
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sujned  that  saturation  effects  on  the  FIR  field  may  be  neg¬ 
lected,  i.e.  S/P  <<  1.  5.  Onlv  homogeneous  broadening  is 

considered.  6.  We  have  also  neglected  r°32/r°i3. 

These  equations  represent  the  response  of  a  three  level 
system  to  an  arbitrary  pump  pulse,  subiect  to  the  assump¬ 
tions  "lust  noted.  An  extremely  simple  form  results  for 
those  conditions  in  which  detuning  on  both  fields  is  absent. 
This  will  be  discussed  in  the  following  section. 

4.  ASYMPTOTIC  LIMITS  AND  COMPARISON  WITH  EXPERIMENT 

4.1  ASYMPTOTIC  LIMITS 

We  first  consider  the  behavior  of  the  pump  pulse  in  terms 
of  a  more  convenient  parameter,  the  absorption,  which  is  de¬ 
fined  as  follows. 

A(z.t)  -  1  -  (45) 

and  a  dimensionless  absorption  length, 

B  *  or°^2  (46) 

We  mav  now  write  egiiation  37  as  follows. 

6  =  ♦  A  -  ln(l  -  A)  f47) 

0 

Thus,  the  absorotion  of  the  pumping  pulse  mav  be  expressed 
as  A(  B  *  "Phis  dependence  is  shown  in  figure  6.  It  is 
to  be  noted  that  the  absorption  is  only  a  function  of  the 


initial  value  of  the  pump  pulse,  which  is  in  itself  a  func¬ 
tion  of  the  retarded  time,  and  the  dimensionless  absorption 
lenqth  6  .  We  will  see  below  that  this  simple  result  com¬ 
pares  well  with  experiment.  Consider  now  the  FIR  field  in¬ 
tensity  as  given  bv  equation  42.  The  algebraic  complexity 
of  this  result  indicates  that  it  is  perhaps  of  value  to  ex¬ 
amine  the  case  in  which  there  is  zero  detuning  i.e.,  *g  = 

=  0.  Thus,  ^12  =  *^2  3“  '^12“’'  »  0®2ir,  and  our 
result  for  the  solution  of  the  FIR  field  intensitv  is  given 


as  follows. 


P  +  1 
o 

P  +  1 


(48) 


where 


5  = 


2 

“32“s 


2 

VI  (t) 

p 


(49) 


We  note  from  ecruation  48  that  for  Pq  >>  1,  If  F  becomes 

small,  that  Ig/IgQ  -  (Po(t  )  )^  •  A  typical  value  of  t  would 

be  of  the  order  of  ten,  so  that  for  these  conditions,  there 
is  a  considerable  steepening  or  compression  of  the  FIR 
pulse. 

We  may  also  define  a  reduced  FIR  gain  coefficient  in 
terms  of  pumo  absorption  A,  the  initial  pump  pulse,  and 

C  • 
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and  this  is  presented  in  figure  6,  Thus,  we  have  an  ex¬ 
tremely  simple  result  for  the  growth  of  the  FIR  intensity 
that  parametrically  depends  upon  the  amount  of  pump  power 
absorbed.  If  we  then  consider  the  situation  of  strong  pump¬ 
ing,  P  >>  1,  then  we  find 


'  (51) 


For  A  =1,  our  results  are  invalid,  since  at  this  point  the 
assumption  of  weak  saturation  on  the  FIR  transition  will 
have  been  violated. 

We  note,  however,  that  not  only  B  but  (or  *q)  and  A 
are  dependent  upon  pressure.  For  the  case  of  the  homogene¬ 
ously  broadened  pumping  pulse  considered  in  the  present 

analysis,  ° the  absorption  coefficient  of  equation  46 
is  proportional  to  pressure.  It  is  therefore  conyenient  to 
define  a  new  parameter  e  that  Is  independent  of  pressure. 


but  proportional  to  pump  input  power  and  cavity  length. 
With  E  =  vfe  may  write  eqation  47  as  A  A(  e  ,  B  ) ,  and 

curves  of  constant  e  correspond  to  curves  of  fixed  cavity 
length,  and  a  fixed  pumping  power.  ’’’he  variation  of  pump 
absorption  with  pressure  occurs  only  through  the  pressure 
dependence  of  the  dimensionless  absorption  length  B  .  See 
figure  7,  we  may  also  use  equations  47  and  50  to  obtain  a 
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This  relation  is  presented  in  figure  8.  Curves  of  constant 

e  are  independent  of  pressure,  and  the  reduced  FIR  gain 
coefficient  is  conveniently  given  in  terms  of  ^  and  the  ab¬ 
sorption  A.  We  see  from  figure  8  that  at  low  values  of  , 
corresponding  to  low  input  power  or  short  cavities,  there  is 
a  broad  maximum  in  the  FIR  gain  coefficient  with  respect  to 
absorption.  This  maximum  shifts  toward  higher  values  of  the 
absorption  as  cavity  length  or  input  power  are  increased. 
We  should  also  note,  that  for  very  high  values  of  absorp¬ 
tion,  our  results  can  not  be  expected  to  provide  quantita¬ 
tive  agreement  with  experiment,  since  the  restriction  that 
S/P  <<  1  will  have  been  violated. 

4.2  COMPARISON  WITH  EXPERIMENT 

In  order  to  obtain  a  comparison  with  experiment  it  is 
necessary  to  obtain  values  for  and  ♦  that  correspond 
to  a  specific  opticaly  pumped  system.  We  consider  the  496 

vm  FIR  CHjf  laser  pumped  by  the  P{20)  9.55  m  CO2 
line. 24-29  Saturation  absorption  of  this  line  has  been  meas¬ 
ured  by  Temkin,  et.  al.  30  ^^d  the  FIR  gain  by  Brown,  et. 
al.31  The  CO2  pumping  pulse  in  the  experiments  have  a  typi¬ 
cal  length  of  100  nanoseconds (FWHM)  and  the  relaxation  times 
are  approximately  10  nanoseconds/p,  where  p  is  the  pressure 
in  torr.  These  conditions  are  consistent  with  the  assump¬ 
tions  made  in  obtaining  our  solutions.  We  shall  now  assume 

that  Tj^  =T2  =  1/  *'^11  “  1/  "'^22*  calculate  the  value  of 


18 


the  absorption  cross  section  for  the  pumping  pulse,  we  note 


where  characterizes  the  homogeneous  linewidth  of  the 

pumping  pulse  ^00^-2  GHz,  and 


Y 
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4-26  X  10^  sec.  ^ 


(54) 


The  remaining  quantity  to  be  determined  in  order  to  obtain 
an  absorption  cross  section,  is  an  average  dipole  moment  for 

for  the  transition  of  interest.  Since  our  simplified 
model  encompasses  only  three  levels,  it  is  necessary  to  sum 
over  the  M  and  K  quantum  levels  in  an  appropriate  fashion. 
This  calculation  is  outlined  in  Appendix  III,  and  we  obtain 
a  value  of  <  >  =  .00213  D,  to  be  substituted  into  equa¬ 
tion  53.  With  this  result,  the  average  of  the  cross  section 
is  given  as  a  =  1.67  x  lO-^^  cra“2  The  initial  state  densi¬ 
ty  of  those  particles  interacting  with  the  pump  pulse  is 
calculated  to  be  =  3.22  x  lO^^  p  cm“3»  where  again,  p 
is  the  pressure  in  torr.  Thus  the  total  absorption  cross 
section  is  given  by  =  5.38  x  10“3  p  cm“^.  These  values 
for  the  state  density  are  taken  at  a  temperature  of  300  K. 
The  dimensionless  intensity,  or  initial  pulse  area  is  given 
as  follows. 
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I 
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However,  the  total  decay  rates  of  level  1  and  level  3  are 
the  summations  of  contributions  from  all  active  sublevels 
weighted  by  the  Boltzmann  factor,  thus,  we  must  sum 

where  fjj^  is  the  statistical  weight  of  the  sublevels.  We 
calculate  an  effective  Boltzmann  constant  f  -  .1118.  Thus 

for  =  10  X  10“^  sec/p,  Y  «  1.118  x  10^  sec/p  and  for 
the  parameters  calculated  above,  we  find  that  ♦q  =  1.44  x 

IO'^Iq/p,  where  Iq  is  the  intensity  in  watts/cm^,  and  p  the 
pressure  in  torr.  From  equation  49,  with  =  .05  D, 

^22  '  1  D,  we  find  the  value  of  C  to  be  24. 

With  a  pump  pulse  of  .5  J  and  a  FWHM  of  100  nanoseconds, 
we  see,  using  equations  46  and  47  ,  a  comparison  of  experi¬ 
ment  and  theory.  This  is  presented  in  figures  9  and  10 
along  with  the  numerical  results  calculated  by  Picharauthu 
and  Sinha.6  Thus,  equation  47  seems  an  adequate  prediction 
of  the  behavior  of  optical  absorption  by  a  pumping  pulse, 
and  only  two  parameters  need  be  obtained  to  produce  a  satis¬ 
factory  agreement  with  experiment.  The  curves  in  figures  9 
and  10  were  obtained  by  utilizing  the  fact  that  °  is  pro¬ 
portional  to  pressure,  and  that  is  inversely  proportional 
to  pressure.  From  equation  49,  we  calculate  the  value  of  ^ 
to  be  ^4  .  Further,  once  having  obtained  the  value  of 

the  absorption  and  6((>^  ,  we  may  compare  the  reduced  gain 
coefficient,  equation  52,  with  the  experiments  of  Brown,  et. 
al.31.  A  slightly  different  value  of  the  absorption  coeffi- 


a 


4 


4 
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cient  was  used  for  this  comparison,  6r&S"o.^2  =  .003  p 

cm“^  and  ®o  “  3*5Po.  It  is  seen  that  agreement  between  the 
theory  presented  here  and  the  experimental  results  for  FIR 
gain  is  satisfactory.  The  close  correlation  with  the  exper¬ 
iment  for  A  ■<  1  perhaps  indicates  that  for  these  low  values 
of  the  FIR  reduced  gain  coefficient,  that  saturation  effects 
do  not  play  a  dominant  role. 


5 .  SUMMARY 

In  the  limit  in  which  the  pumping  pulse  is  long  compared 
with  the  longitudinal  and  transverse  relaxation  times  of  an 
optically  pumped  system,  we  have  obtained  time  dependent  so¬ 
lutions  for  the  absorption  of  an  optical  pumping  pulse  and 
the  subsequent  generation  of  an  FIR  pulse.  The  solution  for 
the  absorption  is  particularly  simple,  and  seems  to  compare 
well  with  several  sets  of  experiments.  The  absorption  de¬ 
pends  only  upon  an  effective  absorption  coefficient,  and,  in 
a  highly  nonlinear  manner,  upon  the  initial  shape  of  the 
pumping  pulse.  Since  we  have  neglected  saturation  effects 
on  the  FIR  wavelength,  the  pumping  pulse  is  not  influenced 
by  the  presence  of  the  FIR  pulse.  This  can  only  be  true  for 
certain  values  of  the  absorption  or  certain  lengths  of  am¬ 
plifier.  The  solution  for  the  pumping  pulse  permits  a  solu¬ 
tion  of  the  response  of  the  FIR  system,  and  by  recasting  the 
equations  for  FIR  field  intensity,  phase,  and  population 
density  differences  in  a  particularly  convenient  form,  it 


has  been  possible  to  solve  this  system  of  equations  paramet¬ 
rically  in  terms  of  the  pumping  pulse  solution.  The  FIR 
field  solution  depends  upon  a  host  of  molecular  parameters, 
upon  the  behavior  of  the  input  pumping  pulse,  and  upon  the 
pumping  absorption  coefficient.  It  also  depends  in  a  com¬ 
plicated  manner,  upon  the  detuning  parameters  of  the  system 
for  both  the  pump  and  the  FIR  signal  fields.  In  the  limit 
in  which  detuning  may  be  neglected,  particularly  simple  re¬ 
sults  may  be  obtained  for  this  complicated  system,  and  these 
results  have  been  shown  to  compare  favorably  with  some  ex¬ 
periments.  The  simple  analytical  form  shows  the  specific 
role  that  the  initial  pumping  pulse  can  have  in  the  shaping 
of  the  FIR  field  intensity,  and  illustrates  how  this  depen¬ 
dence  can  lead  to  pulse  narrowing  effects  as  the  FIR  field 
intensity  grows. 
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Appendix  II 


Appendix  III 


To  calculate  a  value  of  that  we  use  in  our  three  level  model,  we 

must  appropriately  average  over  all  of  the  M  and  K  sublevels  -of  the  tran¬ 
sitions  that  are  capable  of  interacting  with  the  pumping  radiation*  The  single 
narrow  band  absorption  coefficient  of  each  JKll  sublevel  is  given  as 


JKM 


“(e^hcJ^JKM  J 


v(AVj^/2) 


(AVj^/2)‘ 


♦  (v  -  Q(J.K))^ 


(r®  ) 

'■  13''JKM 


where  Q(J,K)  is  the  frequency  of  the  Q  branch  transition  (v^  =  1,  J  =  12,  K 
(Vj  =  0,  J  =  12,  k) ,  and  the  individual  dipole  of  each  JKM  sublevel  is  given  by 


2  hW 
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Freund,  et.  al.  has  determined  y  in  the  above  equation  to  be  0,04686  D. 

The  difference  density  is  related  to  the  total  vibrational 

10  JKM 

state  density  difference  through  the  Boltzmann  factor. 


fr  1 
^^13  JKM 

(r°3) 


^JKM 


and,  using  the  values  for  CH^F  from  reference  33,  we  calculate  at  300  K  that 
f  «  0.1118, 

Since  each  absorption  coefficient  frequency  dependent  through 

the  Lorentzian  lineshape,  the  total  absorption  will  depend  on  the  pump  in¬ 
tensity  at  that  frequency.  An  average  absorption  coefficient 
be  given  by 


We  assume  the  CO2  TEA  laser  has  a  Gaussian  spectral  composition  centered  at 

P(20)»  with  a  bandwidth  (FWHM)  of  Av  =  2GH  . 

C02  ^ 


where  is  the  total  intensity  integrated  over  the  spectral  distribution 

of  the  exciting  CO^  laser  beam. 

The  homogeneous  linewidth  of  each  individual  K  sublevel  is  approximately 
=  40  p  Mlz,  which  is  small  compared  with  the  carbon  dioxide  linewidth  of 
2  gigahertz.  Therefore,  each  individual  K  sublevel  may  be  considered  to  be 
excited  by  a  broadband  pumping  source.  In  the  broadband  pumping  limit  con¬ 
sidered  here,  the  Lorentzian  appearing  in  the  single  narrow  band  absorption 
coefficient  behaves  as  a  6  function  in  the  above  integral  over  the  Gaussian 
frequency  distribution  of  the  pumping  pulse  intensity.  The  indicated  integra¬ 
tion  may  then  be  carried  out.  This  permits  the  identification  of  the  equiva¬ 
lent  level  constants  that  must  appear  in  equation  53. 


2  _  mI  p  f  (2J  ^  1) 


(QP.K)  -  P(20)) 


2 


When  the  values  for  J  =  12,  y  *  .0469  0  are  substituted  into  the  above 

expression,  we  find  that 


=  0.00213  D. 


l  ig.  3.  Schematic  of  three  level  system 


sionlcss  absorption  length  a:  and  dimensionless 


KcJiiccJ  I- 1 1?  (lain  as  function  of  I’orcoiit  Absorption  and 
dimensionless  jx)wcr  parameter  c. 
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Comparison  between  experiment,  numerical  solutions^  and  present  theory,  liq.  47. 
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